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Abstract 

We have collected a 43 event sample of the decay Kl — > e + e~fi + fi~ with 

negligible backgrounds and measured its branching ratio to be (2.62 ± 0.40 ± 

0.17) x 10 -9 . We see no evidence for CP violation in this decay. In addition, we 

set the 90% confidence upper limit on the combined branching ratios for the 

lepton flavor violating decays Kl — > e ± e ± / u T / Li T at B(Kl — ► e ± e ± /r F /r F ) < 

1.23 x 10~ 10 , assuming a uniform phase space distribution. 
PACS numbers: 13.85.Rm, 13.25.Es, 14.40.Aq, 14.80.Ly 
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We present the observation of the decay Kl — > e + e~ fi~ and a measurement of its 
branching ratio. This decay proceeds entirely through the ^7*7* vertex and provides the 
best opportunity for its study. Knowledge of the ^7*7* vertex is crucial in order to extract 
short distance information, including the CKM matrix element Vtd, from Kl — > fi + ^~ 
decays. We are also able to search for CP violating effects in Kl — > e + e~ well 
as the lepton- flavor violating decay Kl — > e ± e ± /i =F / u =F . The previous E799-I experiment 
measured B(K L — > e + e~[i + fi~) = (2. Q^a) x 10~ 9 Q with one event. 

There are several predictions for the branching ratio. Quigg and Jackson |2| have used 
a Vector Meson Dominance (VMD) model to predict a value of 2.37 x 10~ 9 . A phase-space 
model with both CP conserving and CP violating form factors by Uy || predicts values from 
(1.63 ± 0.07) x 10~ 9 for a totally CP conserving decay, to (3.67 ± 0.15) x 10~ 6 for a totally 
CP violating decay. Note that the CP violating form factors increase the branching ratio by 
over three orders of magnitude. This calculation does not take into account any momentum 
dependence of the form factors. An 0(p 6 ) Chiral Perturbation Theory calculation by Zhang 
and Goity § predicts (1.30 ± 0.15) x 10~ 9 . 

The measurement presented here was performed as part of the KTeV experiment, which 
has been described elsewhere ||. The data used were collected during the 1997 run. The 
KTeV experiment, as configured for rare decay searches (E799-II), used two nearly parallel 
kaon beams created by 800 GeV/c protons incident on a BeO target. The kaon decays used 
in our studies were collected in a decay region approximately 65 m long, situated 94 m from 
the production target. 

Charged particles were detected by four drift chambers, each consisting of one horizontal 
and one vertical pair of planes, with typical resolution of 100 fim per plane. Two drift 
chambers were situated on either side of an analysis magnet which imparted approximately 
205 MeV/c of transverse momentum to the charged tracks. The drift chambers were followed 
by a trigger hodoscope bank and a 1.9 m x 1.9 m calorimeter composed of 3100 blocks of 
pure Csl. The fiducial volume was surrounded by a photon veto system used to reject events 
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in which photons missed the calorimeter. The calorimeter was followed by a muon filter 
composed of a 10 cm thick lead wall and three steel walls totalling 511 cm. The first plane 
of scintillators used to identify muons (MU2) was located after 400 cm of steel, behind the 
second steel wall. Two additional 3 m x 3 m scintillator planes (MU3Y and MU3X), located 
after the third steel wall and consisting of one horizontal and one vertical plane, defined the 
acceptance for muons. All muon scintillator planes had 15 cm segmentation. 

The trigger for the signal events required hits in the upstream drift chambers consistent 
with at least two tracks, as well as two hits in the trigger hodoscopes. The calorimeter 
was required to have at least one cluster with over 1 GeV of energy, deposited within a 
20 ns window relative to the event trigger. The muon counters MU3X and MU3Y were 
each required to have at least two hits. In addition, a preliminary online identification 
of K L — > e + e _ /i + /i~ decays required a minimum of three tracks originating from a loosely 
defined vertex. A separate trigger was used to collect K L — > 7r + 7r~7r° decays with subsequent 
Dalitz decays n° — > e + e~7 (K L — > 7r + 7r~7r£,) which were used for normalization. This trigger 
was similar to the signal trigger but had no requirements on hits in the muon hodoscopes 
or clusters in the calorimeter. The preliminary online identification was performed on the 
normalization sample as well. The normalization mode trigger was prescaled by a factor of 
500:1. 

The offline analysis required four tracks from a single vertex. A cut on the vertex 
reconstruction y 2 ensured the four tracks originated from the same vertex. Upstream and 
downstream track segments were allowed at most a 2 mm offset at the bend plane of the 
analysis magnet. The four track decay vertex projected to the calorimeter was required to 
be within one of the beam regions. A track was identified as e ± if it pointed to a cluster 
in the calorimeter with \E/P — 1| < 0.05, otherwise it was identified as a muon. The 
electromagnetic energy resolution of the calorimeter, a(E)/E = 0.45% © 2.0%/ \Je (GeV), 
was determined using a large sample of e ± from K L — > -n T e ± u decays. 

To remove muons which range out in the steel, events in which a muon track had 
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P < 10 GeV/c were discarded. To remove misidentified pions, events in which a muon track 
deposited over 3 GeV in the calorimeter were discarded. Events with excessive energy in the 
veto counters beyond that expected from accidental activity, or extra calorimeter clusters 
not associated with an electron or muon track, were also discarded. 

A large component of the background consists of Kl — > 7r + 7r~7r£> decays in which the 
charged pions simulate muons by punching through the muon filter or decaying in flight. 
A simulation of this background is compared to the data in figure |l|. The invariant mass 
(M eewi ) distributions for the data and for the simulated background are shown after success- 
ful electron identification, but before cuts on vertex quality or extra calorimeter activity. 

The events at high invariant mass are due to two kaons which decayed within the same 
Tevatron RF bucket (double decay events). An event with a double decay contains two 
separate two-track vertices, which usually form a poor four-track vertex. The vertex quality 
cut eliminated nearly all of the double decay events. The cuts on photon veto activity and 
extra calorimeter clusters removed 98% of the background from Kl — > 7r + 7r~7r£,. Figure [IJ 
also shows the data after the vertex quality and extra calorimeter clusters cut. 

Another potentially large background comes from Kl — > /i + /i~7 events in which the 
gamma converts to an e + e~ pair in the vacuum window. A cut requiring that either the two- 
electron invariant mass (M ee ) be greater than 3 MeV/ c 2 or the separation for the two electron 
tracks at the first drift chamber be greater than 3 mm reduced the window conversion events 
by 98.8% while losing only 13.6% of the remaining signal events. We estimate a 0.19 event 
background from the vacuum window conversions. 

The double decay background remaining after all cuts was estimated by examining 
events in which the two electrons were of the same charge sign. Two electrons (or two 
muons) that do not arise from the same decay have the same sign as often as opposite signs. 
Figure |2| shows scatter plots of M eeMM vs. P 2 , where P 2 is the square of the transverse 
momentum with respect to the trajectory from the production target to the decay vertex. 
Figure ||(a) contains signal events with all but M ee ^ and P 2 cuts, while figure ||(b) shows 
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the distribution for events that pass the same cuts, except that the two electrons (and 
muons) are required to have the same sign. There are 30 such like-sign events in the region 
0.25 < M eew < 0.75 GeV/c 2 and P 2 < 0.02( GeV/c) 2 . We assume that this background 
is evenly distributed, and estimate a 0.02 event double decay background within the signal 
region. The M eefJifl vs. P 2 data distribution shown in figure ^|(a) contains 43 events within 
the signal box given by 0.48 < M ee ^ < 0.51 GeV/c 2 and P 2 < 0.00025 (GeV/c) 2 . 

The remaining background at masses below Mk is mainly from Kl — > 7t + -k~7t% decays 
with charged pion decay or punchthrough. Figure |3| shows the fit of the data with all but the 
invariant mass cut to a scaled background simulation, from which we estimate a background 
of 0.03 events in the signal region. The total background in the signal region, from window 
conversions of Kl — ■> /i + /i~7 decays, Kl — ► 7i + 7i^ir t ^ decays in which the charged pions 
either decay in flight or punchthrough the filter steel, and double decay events is estimated 
to be 0.24 events. 

The geometric acceptance of the detector for the signal mode, calculated using a Monte 
Carlo generator which employed the matrix element formulated by Uy |§, was 6.1%. Only 
CP conserving elements of the matrix element were used for the acceptance calculation. 

In order to calculate the branching ratio, we use the Kl — > 7r + 7r~7r^ decay for normal- 
ization. These events were selected by the minimum-bias, two-track trigger described above. 
The cuts applied to the normalization mode were as similar to the signal mode as possible, 
with a few separate cuts specific to this mode, including the requirement that the photon 
cluster in the calorimeter deposit at least 5 GeV of energy in the calorimeter and be > 3 cm 
away from the beam holes. We also required that the invariant mass of the e + e~7 be within 
10 MeV/c 2 of the 7r° mass, and its energy be between 15 and 85 GeV. We determined that 
(2.70 ± 0.08) x 10 11 Kl within an energy range of 20 to 220 GeV decayed between 90 and 
160 meters from the target. The error in this value is dominated by the uncertainty in the 
measured branching ratios for Kl — > 7r + 7r~7r° and 7r° — > e + e~7 ||. 

While similar in most respects to the signal trigger, the normalization mode trigger 
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lacked any muon requirements. One source of systematic uncertainty in the normalization 
involved the response of the muon steel and counters. We evaluated the hit efficiencies of 
the muon counters using high statistics samples of ^ taken in runs with special absorbers 
and magnet configurations. The muon efficiency was ~99% and nearly uniform across 
the counter planes. This efficiency was measured to < 0.5% of itself. We used similar 
runs to study our scattering simulation. By selecting tracks which traversed an overlap of 
adjacent muon counters in MU2, we were able to determine the track position at the counter 
plane independent of the tracking system, and ensure the accuracy of the muon scattering 
simulation for incident muon momenta well below the 10 GeV/c cutoff. 

We estimated the systematic uncertainty in the branching ratio due to any disagree- 
ments between the simulation and data by varying the selection cuts. The contribution 
due to uncertainties in the selection criteria for both the signal and normalization modes 
is 3%. Accidental activity in the detector caused a 3.7% uncertainty in the acceptance, 
while systematic uncertainty due to muon scattering in the lead and steel downstream of 
the calorimeter was estimated to be 0.5%. The change in signal acceptance due to possible 
form-factors in the matrix element of the decay contributed 3%, while the uncertainty due 
to background events remaining in the signal region was estimated at 0.6%. 

The total systematic uncertainty was 6.4%, to be compared to 15.2% statistical uncer- 
tainty. Our final result is thus B{K L -> e + e-//+/i~) = (2-62 ± 0.40 ± 0.17) x 10" 9 , where 
the first error is statistical and the second is systematic. 

Events with like-sign leptons are a signature of lepton flavor violation, and given that 
we found no like-signed lepton events in our sample, we put a limit on the process Kl — > 
e ± e ± /i =F yU =F . Using a flat phase-space generator, similar analysis requirements and signal 
region as that used for the e + e~fi + fi~ analysis, we found that the geometric acceptance for 
like-sign lepton events was 7.3%. With zero events in the signal region, we find that the 
combined B(K L -> e ± e ± fi T fi T ) < 1.23 x 10~ 10 at the 90% confidence level. 

Uy has carried out a phase space calculation for Kl — > e + e~fi + fi~, including a 
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CP violating term in the decay. The form factors used had no q 2 dependence and were 
labelled g<i for the CP even and hi for the CP odd component. According to this model, the 
total measured branching ratio is sensitive to the CP violating term. We measure the ratio 
(fr/hz) 2 < 2.71 x 10" 4 with 90% confidence. 

The angular distributions of the decay products are also a sensitive probe of CP violating 
effects. Define the angle <fi to be the angle between the plane containing the e + e~ and that 
containing the An asymmetry about zero in the distribution sin 20 is unambiguous 

evidence for CP violation. We have observed such an asymmetry in the decay Kl — > 
7r + 7r~e + e~ 0. For the decay Kl — ► e + e~fi + {i~, 18 (25) events are in the positive (negative) 
sin 2(p direction, consistent with zero asymmetry. 

In principle, the K L — ► e + e~/i + /i~ decay is the best mode to investigate the possible 
i^7*7* form factors because there are no exchange terms to complicate the theoretical un- 
derstanding. However, the low statistics makes determination of the form factors difficult 
at present. Figure |] shows the distributions for M ee and M w for the 43 signal events, along 
with the shape expected from a simulation without any momentum dependence in the form 
factors. Although there may be a discrepancy between data and MC prediction in the M w 
distribution indicative of the presence of a form factor, the lack of statistics makes a firm 
conclusion difficult at this time. 

We have measured the branching ratio for the decay Kl — ► e + e~fi + fi~ to be (2.62 ± 
0.40 (stat) ±0.17 (sys) ) x 10" 9 . Note that the VMD model of Quigg and Jackson [§ which 
predicts 2.37 x 10~ 9 is in better agreement with our result than more recent predictions. We 
have placed an upper limit on the CP violating term of Uy JJ, of (f|) 2 < 2.71 x 10~ 4 at 
the 90% C.L, using the branching ratio. We have also searched for lepton flavor violation in 
like-sign lepton decays, and placed a limit B(Kl — > e ± e ± /i =F / u T ) < 1.23 x 10~ 10 at the 90% 
confidence level. The 1999 run of KTeV will soon yield a larger statistics sample of Kl — > 
e + e~fi + fi~ decays. Studies of form factors in these decays should advance the understanding 
of the i^7*7* vertex, thereby constraining one of the largest sources of uncertainty in the 
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determination of the CKM matrix element Vu from Kl — ■> /x/x decays. In addition, a CP 
violating contribution might be detected with a larger sample of K L — > e + e~/x + /i~ decays. 

We gratefully acknowledge the support and effort of the Fermilab staff and the technical 
staffs of the participating institutions for their vital contributions. This work was supported 
in part by the U.S. Department of Energy, the National Science Foundation and the Ministry 
of Education and Science of Japan. 
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FIG. 1. Distribution of M ee/i/x after finding a four-track vertex and identifying two of the 
tracks as for data (hollow circles), and for Monte Carlo simulation of Kl — > 7r + 7r _ 7T£) with 
charged pion punchthrough/decay at the same stage of analysis (dotted line). The simulation is 
normalized to the data below 0.32 GeV/c 2 . The data distribution after cuts on vertex quality 
and any extra clusters in the calorimeter is shown by filled circles. A signal peak is visible at the 
kaon mass. At this stage of the analysis, ~20% of the events in that peak are Kl — > with 
7 — > e + e~ in the vacuum window. 
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FIG. 2. Top plot (a) shows vs. M eeM(U for Kl — > e + e _ /x + // _ , signal events, with all cuts. 
There are 43 events in the signal box. Bottom plot (b) shows same sign lepton events, e ± e ± n T ^ T ■ 
There are no same sign events in the signal box. Note that the area in plot (b) is a factor of 10 
larger than that in (a). 
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FIG. 3. M eewt for data (dots) and the scaled background simulation (line), with Pj 2 < 0.00025 
(GeV/c) 2 . The peak resolution for the signal is approximately 4 MeV/c 2 . 
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FIG. 4. The distributions of M ee and for the 43 events observed in the data (dots) as well 
as that expected from a model without any momentum dependent form factors (line). 
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